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1 Introduction 



Recent remarkable results in iV = 2 extended supersymmetric Yang- Mills theories [1, 2] (for 
further references see these papers) include exact formulas for masses and a theoretical ex- 
planation of basic phenomena such as confinement and chiral symmetry breaking. Obviously, 
it would be very interesting to see how these features emerge in numerical simulations. This 
could contribute to a better understanding of non-perturbative properties of four- dimensional 
supersymmetric quantum gauge field theories (for earlier reviews see [3, 4]). 

A basic property of lattice regularization is that supersymmetry is broken and has to be 
recovered in the continuum limit. In this respect it is quite similar to another basic symme- 
try, namely chiral gauge symmetry [5]. In fact, the authors of ref. [6] showed that in iV = 1 
supersymmetric Yang- Mills theories in the continuum limit there is an intimate connection be- 
tween the restoration of supersymmetry and chiral symmetry. Following their suggestion, in the 
present letter I shall give a prescription on how to perform a supersymmetric continuum limit 
in a broader, non-supersymmetric, renormalizable quantum field theory. As an example, the 
N = 2 supersymmetric SU(2) Yang- Mills theory will be considered, but a similar procedure can 
also be applied in other interesting theories with N = 2 and iV = 4 extended supersymmetry. 

The difficulty of considering supersymmetric theories on the lattice lies in the fact that 
many bare parameters are needed and the parameter tuning is a non-trivial task. In the next 
section the lattice action will first be defined. Then constraints implied by the convergence 
of the path integral will be identified and different minima of the classical potential will be 
classified. Section 3 is devoted to the question of how to tune the bare parameters to reach the 
supersymmetric fixed point in the continuum limit. Another technical difficulty comes from 
the doubling of fermions in popular numerical simulation algorithms. This will be discussed in 
the last section. 



2 Lattice action 

The fields in the iV = 2 supersymmetric SU(2) Yang-Mills theory are: A^^, ipl, Al, B^. 
The gauge field is A^^, /i G {1,2,3,4}, a G {1,2,3}, which is represented on the lattice 
by the SU(2) matrix on links Uxfi = ex.T^{igTaA^^), with the SU(2) generators = r"/2. 
Note that for simplicity the lattice spacing (usually denoted by a) is set to 1 throughout 
this paper, in other words every dimensional quantity is measured in lattice units. The four- 
component fermion in the adjoint representation is described by i^j G {1,2,3}. The 
real scalar and pseudoscalar fields in the adjoint representation are denoted by yl^, r G {1, 2, 3} 
and 5*, s G {1,2,3}, respectively. Instead of the component fields, sometimes, the matrices 
= Till)]., ■^^ = Tj-^^, Ax = TrAl., Bx = TsB^ will also be used. 

The lattice action contains three pieces 

S=Sg + Sf + Ss. (1) 
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The standard Wilson action for the gauge field Sg is a sum over the plaquettes 

^. = /3E(l-^TrC/,,) , (2) 

with the bare gauge coupling given by /3 = 4/^'^. For the fermionic part Sf I shall take the 
Wilson formulation containing the irrelevant parameter r, < r < 1, to suppress fermion 
doublers in the continuum limit: 



r _ ±4 _ _ 



Here /i/ is the bare mass parameter, Kf the hopping parameter, and Ga and Gb are the bare 
Yukawa couplings. Sf is formulated here with Dirac fermion fields. Of course, one can rewrite 
it in terms of two Majorana fermion fields ■^(■^'^^ by the relation = {^^^^ + ii^^"^^) / \^ ■ Then 
the kinetic and mass terms will be multiplied by factors | and the Yukawa coupling becomes 
off-diagonal in the Majorana index. The scalar part of the action is 



±4 ±4 



X IJ,= ±1 IJ,= ±1 

+ ^AAAaA^ + )^AABbAIBI + XbBBbB^ + )^ABAB{Aa:Ba:y J" , (4) 

where the bare mass parameters are fiA,B^ the hopping parameters k,a,b^ and Xaaaa^ • • • 5 ^abab 
denote the quartic couplings. Note that in eqs. (3) and (4) the fields are in general normal- 
ization. In a numerical simulation an appropriate definition of the normalization is given by 
fJ'ip = fJ'A = fJ'B = 1. In general, we have the following relations to the bare masses m^, niA and 
ruB-. 

+ Ar = - — , + 8 = — , + 8 = — . (5) 

Zk^ ka k,b 

Note that the above lattice action of the SU(2) adjoint Higgs-Yukawa model contains all 
renormalizable interactions with the given set of fields if local gauge symmetry and parity 
conservation are required. The corresponding continuum Euclidean action is, in a conventional 
normalization, 

S = l d'x{^F;^{x)F;^{x) + ^{D,A^{x)){D,A^{x)) + ^{D,B^{^^^^ 

+W{xh^D^i;'{x) + ^m\A'{x)A'{x) + ^mlB'{x)B'{x) + m^i^\x)i;' {x) 

+ierstW{x)[GAA'{x) + iGB75B'{x)]i^\x) + Xa[A'{x)A'{x)]' + Xb[B'{x)B'{x)]' 

+ \ab]A'-{x)A'-{x)B'{x)B\x) - \ab)[A'-{x)B'-{x)]' } . (6) 

Here Xa, Xb, X[ab]t X(^ab) are appropriate linear combinations of the quartic couplings in (4). 
The other notations are, as usual, 

F;^ = d.Alix) - d.A^ix) + gerstAl{x)Al{x) , D,A^{x) = d,A^{x) + gerstAl{x)A\x) , (7) 
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and similarly for the covariant derivatives of the other fields. 

The theory defined by the action in (6) is iV = 2 supersymmetric if 

tua = TUB = = \a = Xb = , G\ = Gl = 2\ab] = 2\ab) = 9^ ■ (8) 

(For references and a review see [7].) At zero fermion mass (m^ = 0) there is an SU(2)k (x)U(l)yt 
global symmetry. The SU(2)k part can be seen by writing the action in a (left-handed) Weyl- 
fermion basis by the relations 

M^) = ctcL , i^R = €lC . (9) 

The SU(2)k symmetry transformations act on the left-handed components of the fermion and 
antifermion fields: (■^l^'^cl)- The U(l)yt symmetry is the usual singlet axial symmetry acting 
on fermions and augmented by a chiral rotation of the scalar fields: 

,p{xy = e-'"^^^(a;) , ^(a;)' = ^(a;)e-'"^^ , 

A{xy = cos(2a)yl(a;) - sm{2a)B{x) , 5(a;)' = sin(2a)yl(a;) + cos(2a)5(a;) . (10) 
This turns out to be an anomalous classical symmetry of the action in (6) if 

= , m\ = "^B 5 Ga = Gb , 2\a = 2\b = \ab] - \ab) ■ (H) 

The fermion mass term breaks down SU(2)k to its diagonal subgroup U(l):r, which corresponds 
to the fermion number conservation. In the lattice formulation the Wilson term proportional to 
r breaks SU(2)k (/) U(l)yt to U(l):r even at zero bare fermion mass = 0. The rest of SU(2)k 
has to be restored in the massless continuum limit, similarly to non-singlet axial symmetries 
with Wilson fermions. 

In order that the path integral over the scalar fields with the lattice action (l)-(4) be con- 
vergent, the quartic couplings dominating at large fields have to fulfil the following conditions: 

> AND Ab > AND 

^\ab] > max(0,A(^B)) OR 4A^Ab > max[XfAB]^ {\ab] - \AB)f] } • (12) 

The positivity of Xa,b is in conflict with the supersymmetry condition in eq. (8). Therefore 
at finite lattice spacing one has to take Xa,b > and tune them in the continuum limit to 
zero (see next section). Note that in principle one could stabilize the path integral by adding 
irrelevant higher- dimensional terms to the action, e.g. [A^A^y. However, this would only be 
a complication and could not solve the problem because in the continuum limit the instability 
would reappear. 

The phase structure of the lattice model is, of course, the first important question for non- 
perturbative studies. As a hint for possible phases one considers the minima of the classical 
potential 

V{A,B) = ^mlA'-A'- + ^mlB'-B'- + XA{A'-Ay + XB{B'-By + X[AB]A'-A'-B'B'-X^AB){A'-By . 
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This has the following extrema, with the corresponding values of V, 



A' = B' = 0, Vmin = 0; (14) 
^^.0,B'B'.-^, K..„.-^; (15) 

B' = 0, A' A' = , Kn.,. = ; (16) 



A B = , A A = — ^ — , B B = — ^ — 

A[AB] ~ ^"^A^B "^[AB] ~ ^"^A^B 



4(Af^^j-4A^A^) ' ^^'^ 

(A'-^'')^ = A^A^B'B' , 

^r^r ^ ^\^B - |m^(A[^B] " A(^B)) ^^^^ ^ ^^A^ - |m^(A[^B] " \ab)) 
{\ab] — \AB)y — 4A^Ab ' {\ab] — \AB)y — 4A^Ab ' 

y _ -m\m\{\AB] - \ab)) + m\\B + m|A^ 

4[(A[^^]-A(^^))2-4A^Ab] • ^''^ 

The first minimum in (14) corresponds to the symmetric phase. All the others are defining 
phases with the Higgs mechanism and occur for negative scalar bare mass squares m\,m\ < 0. 
The one in (17) breaks the SU(2) gauge symmetry completely and will presumably not play 
any role in supersymmetry. It will not be considered in what follows. The others break the 
SU(2) gauge symmetry to U(l). The last one in (18) gives two degenerate minima with A and 
B either parallel or antiparallel to each other and, of course, assumes that the denominators 
are non- vanishing. Otherwise, in the general case, there is no corresponding minimum at all, 
except if the bare masses satisfy the condition 

\2 



4A^Ab = (A[^B] - A(^B)) , rrij^^JXB = m^yXA ■ (19) 

In this case the values of the two minima in (15) and (16) are degenerate and there is a minimum 
valley connecting them along the line 



(A'B'Y = A'A'B'B' , A'A'JXa + B'B'JXb = ^ = r= ■ (20) 

If one transforms the vacuum expectation values, say, to the third direction by a global gauge 
transformation, then in the [A^,B^) plane this defines an ellipse. For Xa = A^ the ellipse 
becomes a circle. 

The quantum corrections are, in general, changing the positions of the minima and the 
values of the effective potential at the minima. It can, however, be expected that at least for 
weak couplings the phase structure remains qualitatively the same. Particularly interesting is 
the fate of the minimum valley defined by (20) because of the non-renormalization theorems, 
which tell that the ambiguity of the minimum position is not resolved in a supersymmetric 
situation (see [8, 7]). 
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3 Tuning the parameters 



In this section a way of approaching the supersymmetric continuum limit will be proposed. 
Only the Higgs phases will be considered here because they are more interesting physically. 
Since the N = 2 supersymmetric Yang-Mills theory is asymptotically free, the perturbative f3- 
functions can give qualitative information about the behaviour of the lines of constant physics 
(LCPs) in bare parameter space. These lines, which are defined by keeping the independent 
renormalized couplings fixed, go to the continuum limit, if it exists. The naive prescription to 
simply look for the LCP in bare parameter space, which corresponds to the supersymmetric 
relations in (8), does not work because of the conflict with the path integral stability conditions 
(12). 

This means that one has to consider the renormalization group equations (RGEs) of all 
renormalizable couplings, including also those which are not asymptotically free. The two-loop 
RGE for the gauge coupling is: 



dt 



(47r)2 (47r)4 

For the other couplings [9, 10] one can give the equations immediately in terms of the ratios 



(21) 



R 



G 



A,B 



A,B 



A 



rA,B 



A,B 
2~ ' 



A 



-[AB] 



[AB]_ 
2 ' 



A 



r(AB) 



(AB) 



(22) 



Up to one-loop order we have 

dRA ^ J_ 
dt (47r)2 

dRB _ 
IT ~ (47r)2 



• ^Ra{R\ - 1) + C 



^Rb{R\ - 1) + C 



(47r)^ 



(47r)^ 



dvA 
dt 

drs 
dt 



{Airy 



[3 - WvA + WvaRI - 4:R\ + 88rl + 2rlAB) + ^rf^B] " '^nAB)r[AB]] + O 
[3 - + WtbRI - 4:R% + 88rl + 2rlAB) + ^^^b] " ^nAB)r[AB]] +0 



(47r)4 
^4 



dr 



[AB] 



dt 



{AttY 



[3 - 16r[AB] + 8r[AB]{Rl + R'b) " ^R^Rl 



+8{rA + rB){5r[AB] - r(AB)) + ^Qr^[AB] + ^'^UB)] + ^ 



{AttY 



dr 



(AB) 



dt 



{AttY 



[-3 - Wr^AB) + 8r^AB){R\ + Rl) 



+ 16(r^ + rB)r(^AB) - 20r(^5) + 'i2r(^AB)r[AB^ + O 



(47r)^ 



(23) 
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Here the right-hand sides vanish at the supersymmetric point 



susy 



Ra = 1, Rb = 1, ta = 0, rB = 0, r[AB] = ^, r^AB) = ^ 



(24) 



In other words, supersymmetry corresponds to a fixed point of the coupling ratios, which we 
shall call in the following simply the supersymmetric fixed point. The stability properties of 
this fixed point play a crucial role in parameter tuning for the supersymmetric continuum limit. 

The RGEs in eqs. (21) and (23) hold, for instance, for the bare couplings at fixed renor- 
malized couplings. This means that they give the LCPs in bare parameter space. In this 
case t = ^og[m~^y^-^^i) , where rUphysicai is some physical mass in lattice units. The continuum 
limit corresponds to ^ — )• oo. The same one-loop equations also describe the change of the 
renormalized couplings for fixed bare couplings, if the right-hand sides are multiplied by overall 
minus signs. This means that in the latter case the flow of couplings is reversed: what was an 
attractive fixed point before becomes a repulsive one and vice versa. 

The stability of the supersymmetric fixed point can be investigated in the linearized ap- 
proximation of eqs. (23) near r = r.. 



susy ' 



dr 

dt 



(47r)2 



D(r 



susy 



) + 



(25) 



The matrix D is given by 



D 



/ 


16 














\ 







16 
















-16 











4 










-16 








4 







-8 


-8 


16 


16 


16 


4 


v 


8 


8 


8 


8 


16 


-4 / 



(26) 



The eigenvalues ai...6 and eigenvectors ei . e of this matrix are the following: 



a2 = 
as = 
0.4 = 

05 = 
ae = 



= 
16 
16 

-8 
-4 
24 



ei = (0,0,1,-1,0,0) , 

62 = (1,1,0,0,4,4) , 

63 = (2,0,-1,1,4,4) , 

64 = (0,0,1,1,-2,4) , 

65 = (0,0,1,1,-1,-3) , 

66 = (0,0,1,1,6,4) . 



(27) 



Since there are eigenvalues with different signs, the supersymmetric fixed point is neither at- 
tractive nor repulsive: its behaviour depends on the direction of approaching it. Three of the 
eigenvectors, namely 61^3^4, point towards directions where the conditions for path integral sta- 
bility in (12) are violated. They have to be excluded by appropriate parameter tuning. Two 
out of the remaining three directions (62,5) have positive eigenvalues and there is only a single 
direction (65) with negative eigenvalue. A consequence of this is that the supersymmetric con- 
tinuum limit at [g"^ = 0, r = Vsusy) can only be reached by an LCP if it comes from this unique 
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phase structure 
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Figure 1: The schematic view of the generic phase structure in the plane of scalar 
hopping parameters near the critical point n'^ . Besides nA,B all other bare parame- 
ters are fixed; S is the symmetric phase, A and B correspond to non-zero vacuum 
expectation values of the A- and 5-fields, respectively. In the phase AB both scalar 
fields have non-zero vacuum expectation values. 

direction in bare coupling space. All other LCPs go either to infinity in some component of 
r, or end by violating eq. (12). The LCPs reaching a supersymmetric continuum limit are the 

ways to super symmetry. 

The remaining question is the tuning of the bare mass parameters [m^,mA,mB)- In a 
numerical simulation they are usually represented by the corresponding hopping parameters 
(k^, KjI, Kb). For fixed bare couplings [g"^ and r) one expects the existence of a critical line 
k,'^[k,^) = [k,a = k,a{'^iI>)t i^B = K'si'^ip)) where all scalar boson masses are zero in lattice units. 
The existence of such a line of critical points, corresponding to a second-order phase transition, 
has to be checked by numerical simulations. The condition for tuning the fermion hopping 
parameter k,^ is the restoration of the global SU(2)k symmetry. This happens if the diagonal 
part of the fermion mass matrix (in triplet index) vanishes. Denoting this value of k,^ by k^, 
the remaining task is to tune the scalar hopping parameters k,a,b to the critical point k,'^ = 
(k^, k^, k,b)- Going to some other point of the critical line k,'^[k,^), the result is supersymmetry 
softly broken by a Dirac fermion mass. 

On the basis of the analysis of the minima of the classical potential the generic phase 
structure near k,'^ is expected to look as is shown in fig. 1. In sector AB the relevant minimum of 
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the effective potential corresponds to the minimum in eq. (18). As has been discussed there, for 
4Aj4Ab = (A[j1b] — A(j4b))^ this sector is shrinking to a line defined by (19). This refers, however, 
only to the classical theory. In the case of interest for supersymmetry, when Xa = the points 
of the minimum valley defined in (20) can be parametrized by the global chiral transformation 
of the scalar fields given in (10). This symmetry is, however, anomalous; therefore it has to 
be explicitly broken in lattice regularization. In fact, the Wilson term proportional to r in the 
lattice action (3) breaks it and the anomaly is emerging in the continuum limit in the well- 
known way [11]. Since the chiral symmetry is broken in the fermion sector, the degeneracy of 
the minima in the valley is removed by quantum corrections and the AB phase never collapses 
to a line. Therefore the ratio of the vacuum expectation values of the A- and 5-fields can be 
arbitrarily tuned by the angle of approaching the critical point k,'^ in the (k^i, k^)- plane. This 
ratio is a relevant parameter of different supersymmetric continuum limits. Before discussing 
the tuning of the overall lengths of the vacuum expectation values let us return to the RGEs 
in (21), (23). 

As stated before, the RGEs with a minus sign on the right-hand sides can also be used, for 
sufficiently small couplings, to describe the change of the renormalized couplings for fixed bare 
couplings. Since the signs of the eigenvalues in (27) are now reversed, there are two attractive 
directions: 62 and ee. Tuning the initial values of the renormalized couplings with the help 
of the bare couplings to the plane spanned by these two vectors, and then going towards the 
critical point for fixed bare couplings, the ratios of renormalized couplings are approaching the 
supersymmetric relations in (8). 

The overall ratios of the renormalized vacuum expectation values of the A- and 5-fields 
to the A-parameter of the gauge coupling vat I ^latt-, vbt/ ^latt can be tuned as follows. 
Proceeding as described in the previous paragraph and assuming that one starts from a region 
where the scalar masses are of order 1 in lattice units, the initial values of these ratios are 
essentially fixed by the initial value of the bare gauge coupling g"^. Since Alatt exp(— 27r^/5'^) 
becomes smaller for decreasing g"^, the ratios vat I ^latt-, vbt/ ^latt become larger. That is 
for large f3 = ^g~^ one obtains large vacuum expectation values in units of the A-parameter of 
the gauge coupling. Note that, in contrast to pure gauge theory or QCD, in Higgs phases the 
numerical simulation at large (3 is not a problem [12]. 

In summary, the renormalized vacuum expectation values of the scalar fields can be arbi- 
trarily tuned by the two scalar hopping parameters. Therefore the points of the moduli space 
of quantum vacua investigated in ref. [1] correspond to different physical theories, much as, for 
instance, QCD with different values of the quark masses. 

4 Outlook 

Tuning the parameters in the SU(2) adjoint Higgs- Yukawa model on the lattice to reach the 
N = 2 supersymmetric continuum limit can be guided by the fixed point properties of the 
renormalization group equations. In particular, one can combine the use of a single attractive 
direction for bare couplings and two linearly independent attractive directions for renormalized 
couplings. For the successful tuning of bare parameters of the lattice theory to reach super- 
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symmetry in the continuum limit a decisive question is the non-perturbative phase structure 
(see fig. 1). If the appropriate phase structure is reproduced, the points of the quantum mod- 
uli space of vacua are the endpoints of a two-parameter family of LCPs, which represent all 
possible ways to iV = 2 supersymmetry in the SU(2) adjoint Higgs- Yukawa model. 

The lattice methods are well suited for these kinds of investigations. There is extensive 
experience in simulating SU(2) Higgs and Yukawa models with scalar and fermion doublets. 
(For the more difficult case of Yukawa models see e.g. [13].) The replacement of doublets by 
triplets in the adjoint representation does not seem to be a problem. 

The crucial question, which finally will decide whether the numerical simulation can be 
performed by presently known methods, is the behaviour and implementation of the fermionic 
determinant. In several popular and effective simulation algorithms, such as the Hybrid Monte 
Carlo, the number of fermion fields has to be doubled, which would destroy the balance between 
fermionic and bosonic fields. This is, of course, prohibited in a supersymmetric continuum limit. 
There are algorithms, such as the hybrid classical-Langevin algorithm (see e.g. [14]) where the 
appropriate number of fermions can be simulated, but only under the condition that the phase 
of the fermionic determinant is negligible. Taking into account small fluctuations of the phase 
is possible by the method developed in [15]. Large phase fluctuations are, however, at present 
not tolerable. 

In case of the SU(2) adjoint Higgs- Yukawa model the fermion determinant is real therefore 
the phase factor can only be ±1. Previous experience with fermion determinants in other models 
makes it plausible that the fluctuation of this sign is not important. In any case, this can be 
checked by explicitly calculating the fermion determinant [15]. If the effect of the sign of the 
fermion determinant is small indeed, then the numerical simulation of iV = 2 supersymmetry 
becomes possible. I hope to return to this question in a future publication. 
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